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form factor F^^^^* . Numerical results are presented for the intermediate energy kinematics of 
HADES experiments. 
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I. INTRODUCTION 



The ri meson as member of the octet of Goldstone bosons has the valence quark struc- 
ture {uu + dd — ss)/y/3 when choosing the mixing angle of 19.5" in superimposing the 
octet-778 and singlet-Tyo- Various conservations laws forbid low-order decays causing the 
very narrow rj width. This makes the rj decays sensitive for testing invariances of the 
standard model. The hidden strangeness content lets one argue for some sensitivity to 
the strangeness content of the nucleon when considering 7] production off nucleons. Con- 
sequently, the 77 production and various special decay channels were subject of intense 
investigations since some time, both experimentally and theoretically. 

There is a rich data basis for 77 production in nucleon- nucleon collisions providing a test 
ground for meson production in strong interaction processes, in particular near threshold. 
Due to the iso-scalar character of the rj meson, rj production off the nucleon proceeds via 
selected baryon resonances thus allowing differential access to resonance properties. 

Furthermore, the Dalitz decay rj ^e^e~ constitutes a prominent source of di- 
electrons in intermediate-energy heavy-ion collisions. Indeed, the recent HADES data 
[1] exhibit a sizeable yield of e^e" in the invariant-mass region 150 - 500 MeV which is 
essentially attributed [2] to 77 Dahtz decays, but A Dahtz decays and non-resonant vir- 
tual bremsstrahlung [3] contribute in this region, too. The primary aim of the HADES 
experiments [4] is to seek for signals of chiral symmetry restoration in compressed nu- 
clear matter. For such an enterprize one needs good control of the competing background 
processes, among them the mentioned 77 Dalitz decays. 

77 Dalitz decays depend on the pseudo-scalar transition form factor. Such form factors 
encode information on hadrons which is accessible in first-principle QCD calculations or 
abridged variants thereof, such as effective hadron theories or QCD sum rules. In so 
far, experimental information on transition form factors is quite valuable [5]. Given this 
motivation we consider here the process of Dalitz decay of the pseudo-scalar r] meson 

77 ^ 7 7* ^ 7 e~ e"*", (1.1) 

where 7* denotes a virtual photon. Obviously, the probability of emitting a virtual pho- 
ton is governed by the dynamical electromagnetic structure of the dressed transition 
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vertex 77 — > 77* which is condensed in the transition form factor F^-y^* . If the decaying 
hadron were point hke, then a calculation of mass distributions and decay widths would be 
straightforward along the standard quantum electrodynamics (QED). Deviations of the 
measured quantities from QED predictions directly reflect the effects of the form factor 
and thus the internal hadron structure. 

Often, the production process and the decay process are dealt with separately. With 
respect to available new data from HADES [6] the reaction 

P1+P2 ^ Pi + P2 + ^ ^ p'l + P2 + 7 + e"^ + e", (1.2) 

which will be improved in near future, we consider the complete reaction (1.2) in which 

(1.1) figures as a subreaction. The employed framework is that of an effective description 
in a hadronic basis. To be specific, we are going to parameterize the 77 production subre- 
action in nucleon-nucleon (NN) collisions within the one-boson exchange (OBE) model. 
Such an approach has been utilized fairly successfully by various authors in the past. For 
instance, in Ref. [7] a model based on OBE with nucleons only has been facihtated for 
pp —>■ pprj and pd —>■ pdrj reactions. A non-relativistic approach was proposed in Ref. [8] . 
In [9] a description based on OBE without internal meson conversion and with the reso- 
nance S'ii(1535) has been elaborated for A^A^ NNrj in a region sufficiently above the 
threshold. A detailed analysis of NN — > NNrj has been worked out by Nakayama and 
collaborators [10]. For a study of the role of nucleon resonances in rj photo-production we 
refer the interested reader to Ref. [11]. 

In contrast to a factorized description with production of an on-shell rj and an indepen- 
dent decay of an on-shell 77 we attempt here a complete description of the whole reaction 

(1.2) . That means we supplement the subreaction A^A^ NNi] by the Dalitz decay 
part thus dealing with intermediate off-shell 77 meson. This allows considering di-electron 
masses larger than the pole-mass of the 77 meson (for the on-shell 77 meson the invari- 
ant mass of the di-electron is restricted by the mass of the 77 meson). This can furnish 
additional information on the transition form factor in a larger kinematical region. 

Our paper is organized as follows. In section 2 we introduce the rj — » 77* transition 
form factor. Section 3 is devoted to the theoretical background for dealing with the reac- 
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tion pp — > ppri ppqe^e" . It is essentially based on an extension of the effective model 
[12, 13] adjusted to vector (V) meson production in NN — > NNV reactions. The model 
utilizes a direct calculations of the relevant tree-level Feynman diagrams within a phe- 
nomenological meson-nucleon theory. The model parameters have been previously fixed 
from independent experiments and adjusted to achieve a good description the available 
experimental data [12, 13]. However, in the present paper also diagrams with excitation of 
nucleon resonances with masses close to the mass of a nucleon plus 77 meson are included. 
These are jS'ii(1535) , Z)i3(1520) and i^ii(1440) resonances. The corresponding effective 
constants, whenever possible, are obtained from the known decay widths of direct decay 
into T] channel or radiative decay with subsequent use of vector meson dominance. We use 
effective constants commonly adopted in the literature and obtained from different con- 
siderations, e.g., SU(3) symmetry or adjustment to photoabsorbtion etc. [10]. Numerical 
results are presented in section 4, where we consider separately the reaction NN NNrj. 
Results for the full reaction (1.2) are described in section 5 with emphasis on the role of 
the 7777* transition form factor. The conclusions are summarized in section 6, and some 
formal relations are relegated to the appendices. 



II. DALITZ DECAY AND TRANSITION FORM FACTOR 

Let us first consider first the process of a two-photon decay of an 7] meson. The effective 
Lagrangian describing the vertex 77 — > 77 reads [14-17] 



where is the electromagnetic four-potential, denotes the pseudo-scalar 77 meson 
field and is the corresponding coupling constant. The fully antisymmetric Levi- 
Civita symbol e^i^ap is normalized as 60123 = —1- The 77 decay width follows from (2.1) 

as 

s3/2 

and serves for a determination of the coupling constant f^rj-yy from experimental data. The 
square of the 77 invariant mass is denoted by s^, = for an on-shell 77 meson. (Note 



that contrarily to the vector meson case [18], instead of a factor 1/3 [due to averaging over 
three projections of the spin of the vector particle] in eq. (2.2) a factor of 1/2 appears due 
to two photons in the final state.) Experimentally, the branching ratio Br{r] — > 77) = 
r^_^^/r^,toi is known as (39.38 ± 0.26)% [19]. Eq. (2.2) yields |/^^^| ~ 0.025 Ge^-^ for 
the known total width Tn^tot = (1.3 ±0.07) keV. Since in our further calculations the cross 
section is directly proportional to |/^^^ the sign of the coupling constant does not play 
a role and for definiteness it has been taken positively. 

In the decay (1.1), however, one of the emitted photons is virtual with a time- like four- 
momentum and, consequently, the Lagrangian (2.1) must be supplemented by inclusion 
of the corresponding transition form factor (FF). We employ the following procedure: 



where s^* is the di-electron invariant mass squared and frpfjiO) = |/?777l- Formally, eq. (2.3) 
can be considered as the definition of the transition form factor. By a direct calculation 
of the corresponding diagram for the decay rate dV/ds^* for the f] meson one finds [5] 



with a as electromagnetic fine structure constant and //g the electron mass. In the 
kinematical region we are interested in the terms with /i^ can be neglected. Putting 
|F^-y^* (s-y*)|^ = 1 would mean neglecting the finite size and internal structure of 77. It 
is seen that the differential decay width dT^^^e+e- /ds^* is determined by (i) a purely 
kinematical (calculable) factor, (ii) the real photon decay vertex 77 — 77 (known from 
experimental data), and (iii) the (wanted) transition FF F^-y-y* (s^*). Hence, eq. (2.4) evi- 
dences that by measuring the invariant mass distribution one can get direct experimental 
access to the transition FF [5, 20-22]. For a on- mass shell rj meson, the value of the 
di-electron mass Sj* is kinematically restricted by Sj* < m^. In case of reactions of the 
type (1.2) the intermediate r) meson can be off-shell and, in principle, the value of the 
di-electron invariant mass can be larger than the r] pole mass. This situation will be 
investigated below. 

A quite successful theoretical approach to EE's is based on the vector meson dominance 
(VMD) conjecture [14, 17]. Reasonably good descriptions of elastic EE's in the time-hke 




(2.3) 




(2.4) 
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region has been accomplished, indeed. By using the current-field identity [14] 




(2.5) 



with coupling constants f^p and f^^ known [23, 24] from experimentally measured elec- 
tromagnetic decay widths, one can also compute the transition form factor F^^^J^ [s^*) by 
evaluating the corresponding Feynman diagrams (see below). Contrary to the transition 
FF of vector meson production [18], the 77 meson FF, computed within such an approach, 
exhibits a rather good agreement with previous data [5, 20-22]. 

III. MODEL 

We implement the discussed Dalitz decay 7] 7e^e~ into a more general process of 
di-electron production in NN reactions with intermediate 77. Consider the reaction (1.2) 
for which the process (1.1) enters as a subreaction. The invariant cross section is 



where the factor 1/n! accounts for n identical particles in the final state, \T\^ denotes 
the invariant amplitude squared and dr/ is the invariant phase volume which is chosen 
within the so-called "duplication" kinematics [25], i.e. the one which exploits invariant 
two-dimensional phase volumes R2 describing the decay kinematics of a real or virtual 
particle with the invariant mass squared s (s > 0) into two particles, which can also be 
either real or virtual. This kinematics is schematically depicted in Fig. 1. Such a choice 
of kinematical variables is extremely useful if one considers specific classes of Feynman 
diagrams which allow to separate some vertices in a factorized form. Then in the total 
cross section some integrations can be performed analytically (see [3, 18]). For the present 
task we need to consider only such types of diagrams which allows to factorize the meson 
decay vertex from the parts describing the creation of the meson in a NN interaction. 
Consequently, analytical integrations can be executed over the variables connected with 
the decay vertex. 
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A. Nucleon Current 



The invariant amplitude T is evaluated here within a phenomenological meson-nucleon 
theory based on effective interaction Lagrangians which include scalar (cr), pseudo-scalar 
isovector (tt), neutral pseudo-scalar (77) and neutral vector {uj) and vector isovector (p) 
mesons (see [3, 10, 12, 13]) 

CaNN = q^nnNN^^, (3.2) 
C^NN = -^N^5rd,{T^n)N, (3.3) 

Cr^NN = -^-^N^,i^d,%N, (3.4) 

CpNN = -gpNN (N-f^rN^,^ - -^Na^^rNd^^p^^ , (3.5) 

C^NN = -g.NN (iV7^iV^>^ - ^^^M-^^'^^) ' (3-6) 

where N and # denote the nucleon and meson fields, respectively and bold face letters 
stand for isovectors. All couplings with off-mass shell mesons are dressed by monopole 
form factors Fm = (A|^ — / (Af^ — /c|^), where k\j is the four- momentum of a virtual 
meson with mass ^m- 

The Lagrangians (3.2 - 3.6) are needed in evaluations of Feynman diagrams describ- 
ing the Dalitz decay of the 77 meson created from nucleon bremsstrahlung due to A^'A'" 
interaction via one-boson exchange, see Fig. 2a. 



B. Internal Conversion Current 

The 77 meson can also be produced by an internal conversion of the exchanged mesons, 
the so-called conversion current. The dominant exchange vector mesons {V) in this case 
are uj and p mesons with the interaction Lagrangians 

App = (3.8) 

The corresponding diagrams are exhibited in Fig. 2b. 
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The VVt) vertices in the conversion diagrams have been calculated from the radiative 
decay V ^ rj^ within the VMD model [10, 26]. Correspondingly, the vertex form factor 
VVrj is chosen as 

Fyy,{A, kl, kl) = ^2 _ fcf ^2 _ fc2 (3-9) 

which, in accordance with the procedure of determining the coupling constant, is normal- 
ized to unity when one vector meson is on-mass shell and the other one becomes massless, 
e.g., Fvvr,{K kl = ml, k^^O)^ 1. 



C. Nucleon Resonance Current 



In the threshold-near kinematics for rj production in NN reactions there are a few 
nucleon resonances with masses m^r* ~ mi\j+m^ which can contribute to the cross section. 
These are S'ii(1535) and Di3(1520) with odd parity and spins | and |, respectively, and 
the spin- 1 even-parity Pii(1440) Roper resonance with nucleon quantum numbers. The 
corresponding interaction Lagrangians for 5'ii(1535) are [10, 27] 

AiViVl,35 = ^^^^^1535 N*^^Qf^<^^N + h.c, (3.10) 
7T7.JV* — ^AT 

jC.nn,,s. = N*^ d^T^^N + h.c., (3.11) 

>C.^iVi535 = ^^^^^1535 N*^^a,.d''^^^N + h.c., (3.12) 

niN* + TTlN 

^pNN,s3s = ^'"^"I"" 7V*75a^.a-T#^7V + /,.c.. (3.13) 

niN* + TTlN ^ 

For the spin-| resonance Dis{1520) we employ 

AiViv,5.o = ^^'^^^N*^e,A^h5d''%N + h.c., (3.14) 
J^.NN,,,o = ^^^^^N*^e,.izh,d''T^^N + h.c., (3.15) 

C.Nm,,o = ^^^|^A^*^e^.(^)7a^"'^A^ - ^^^^^^do.N*^u;''^N + h.c, (3.16) 

(1) (2) 

= ^^f^iV*^e^.(^)7aTp"'^A^ - ^-lf^daN*'^Tp''''N + h.c., (3.17) 

where 6^1,(2;) = — (A(l + 42;)/2 + 2;)7;,7,^. The field strengths uj°'^ and are ^"^{x) = 
cP^^^{x) — d'^^'^{x) and p"''{x) — 9"#^(a;) — 9''<?^(a;), respectively. The interactions (3.14 



T'f (P) = g>^^ - --i^Y - ^^z^ - Tjz— (7'^^' - l^'P^) , (3-19) 



- 3.17) with Q^v{z) have been extensively discussed in the Uterature (see, e.g., Ref. [28] 
and further references quoted therein). The form of ©^1/(2;) is chosen in such a way that 
the respective interaction Lagrangian obeys the same point transformation as the free 
Lagrangian, which ensures that the matrix elements are independent of the parameter A, 
usually taken as A = —1. The off-shell parameter z remains free thereby. We employ 
here z = —1/2. Also, the choice of the form of higher spin s > | propagators has 
been a subject of discussion in the literature [3, 29-32], e.g., with concerns about the 
spin-projector operator P| (off-mass shell vs. on-shell) or the ordering of the product of 
energy projection operator P^* -|-mjv* and (only on the mass shell these two operators 
commute). In the present paper we take the spin-| propagator in the form [3] 

sr{p) = -i^^^^^vnn (3.18) 

2 1 TJlj^i^ 2 

where the spin projection operator is of the form 

1 ^ ^ 2P^'P'' 1 

as commonly adopted within the Rarita-Schwinger formalism [33]. 

For Pii(1440) the effective Lagrangians are similar to that for 5*11, except for the 
relative signs in 77 and tt couplings and the aNN* Lagrangian, i.e. 

jO,nn,,,o = - ^^"^"y"" N*757,d^%N + h.c., (3.20) 

rriN* + rriM 

JC.nn,,,o = —^^^^^^N*^5l^d^T^^N + h.c, (3.21) 
niN* + niN 

>CaJVJVi44o = gaNN,,,oN*N$„ + h.c., (3.22) 

j0.nn,,,o = N*a,,d''<l>^^N + h.c., (3.23) 

jO,nn,,,o = ^P^^i44o ATV^.a-T^^Ar + fe.c.. (3.24) 

These interaction Lagrangians enter the nucleon resonance current exhibited in Fig. 2a. 

In order to account for the finiteness of the resonance widths F, the resonance mass in 
the corresponding propagators is augmented by an imaginary part, mjv* — > mjv* — irjv*/2. 
Also for the resonance and conversion currents all the vertices with off-shell hadrons are 
dressed by form factors 

A"^ 

F{^,P^) = T2 — r-9 (3-25) 
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D. Effective parameters 



The effective parameters of the Lagrangians (3.2 - 3.6) which determine the pure nu- 
cleoli current diagrams are basically the ones from the OBE Bonn potential [34] which 
have been used in our previous calculations for vector meson {cu, p and 0) production in 
NN interactions [12, 18], see Tab. I. The nucleon cut-off parameter for the bremsstrahlung 
NN-q vertex is Af^iv = 1-2 Gel/. 

TABLE I: Coupling constants and cut-off parameters for the nucleon current [34]- 

9MNN ^MNN [GeV] 

TT /ttATAT = 1.0 1.3 

r/ /r;iVJV = 1.79 1.8 

(T QaNN = 10. 1.8 

P 9pNN = 3.5 1.3 

(Kp = 6.1) 
^ QuiNN = 15.85 1.5 

(^c = 0.0) 

The effective coupling constants for the nucleon resonance currents, whenever possi- 
ble, have been obtained from the known decay widths of direct decay into rj channels or 
radiative decay with subsequent use of VMD. Otherwise, we use effective constants com- 
monly utilized in the literature and obtained from different considerations, e.g., SU(3) 
symmetry, fit of photo- absorbtion reactions etc. (see, e.g. [10]). The few remaining 
less known cut-off parameters are taken either close to the ones from OBE potential 
(for instance, the ApjVAr* and A^j^^* cut off's are chosen equal to A^^jvat = 1.5 GeV; 
A^jvATisas = A^iVAri44o = A^atat = 1.3 GeV, cf. Tab. II) or adjusted to experimental data 
(see below). The resonance cut-off parameters for the bremsstrahlung N*Nr] vertex have 
been taken as A-'^NNi^an = 1-3 GeF, A^^j^^^^^ = 1.1 GeV and A'^^^^^^ = 1.2 GeV, respec- 
tively. 

The coupling constants gft^^^^ and gpprj for the conversion current have been calculated 
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TABLE II: Coupling constants and cut-off parameters for the nucleon resonance current. For 

the spin-^ resonance Di^ the off-shell parameter z = —1/2 is used and the second coupling 
constant, eq. (3.1 7), with vector mesons is given in parenthesis. For a reasoning of g^rNN* and 
gr)NN* for N* = Z)i3(1520) see Appendix A. 



Sii(1535) r>i3(1520) Pii(1470) 

QMNN* ^^[GeV] gMNN' A.[GeV] guNN* A [GeF] 



TT 


1.25 


1.2 


1.55 


1.0 


6.54 


1.3 




2.02 


1.3 


8.3 


1.2 


0.49 


1.3 


U 


-0.72 


1.5 


-2.1(0.7) 


1.5 


-0.37 


1.5 


P 


-0.65 


1.5 


6(-2.1) 


1.5 


-0.57 


1.5 



[10] from a combined analysis of the radiative decay within the VMD model and with a 

SU(3) effective Lagrangian which provides 

g^^r, ^ 4.85; ~ 4.95. (3.26) 

(Note that a naive direct calculations of these constants within the VMD model can 
provide slightly larger values gvvr) ~ 6.0.) The corresponding cut-offs are A^^^^^ = A^^^ = 
1.6 Gey. 



E. Cross section for pp — > pprj 

The invariant amplitude T in eq. (3.1) corresponding to the Feynman diagrams in 
Fig. 2 can be cast in a factorized form 

where the amplitude Tj^\[^NNr) describes the process of production of an off-shell rj meson 
in a A'^A'" collision, while the amplitude T^ J,^g+g- describes the Dalitz decay of the produced 
7] meson into a real photon and a di-clcctron. In the propagator of the rj meson the mass 
has been replaced by — to take into account the finite life time of the r] 

meson. 
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As mentioned, for such factorized Feynman diagrams one can separate, in the cross 
section, the dependence on the variables connected with the Dahtz decay vertex and 
perform the phase space integration analytically. Equation (3.27) allows one to rewrite 
the differential cross section (3.1) in a factorized form as well, 

- ^IP |2r (Sv--h*f I ,5 tot /oooN 

~ 4vr (y^ - m,f + \Vl,^, ^nn^nn,, (3.29) 

where the integral over the final di-electron and photon variables has been performed 
analytically (see for details [3, 18]). dV^^^e+e- / ds.y* is defined by eq. (2.4) but with 
— > and the production cross section of a pseudo-scalar meson with 77 quantum 
numbers but with 7^ is 

" '^NN^NNr) 



X \T''Nl^NNn?dsi2dR2{N^N2^ Sr,s^2)dR2{si2^ N[N'^),{2>m) 



spins 



where the two-particle invariant phase space volume R2 reads 



fl,(„6^ed) = i%^!?^dSi:. (3.31) 

In what follows we are interested in production and subsequent decay of intermediate 77 
mesons, thus being off-mass shell, 7^ m^. We assume that for the off-shell 7] the net 
coupling constant /^/^^(O) in eq. (2.3) obtained from the experimental data via eq. (2.4) 
remains the same. The off-mass shellness of 77 mesons is taken into account by including 
additional effective form factors in the corresponding vertices. 



IV. T) PRODUCTION CROSS SECTION 

It can be seen from (3.29) that the peculiarities of the cross section for the full reaction 
NN NN'^e^e~ are basically determined by the subrcaction NN — >■ NNt] with creation 
of a virtual 77 meson. Hence, before analyzing the full reaction, we proceed with a study 
of the subreaction NN — > A'"A'"77 for production of an on-shell 77 meson. 
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A. Initial and Final State Interactions 

It ought to be mentioned that the Feynman diagrams depicted in Fig. 2 cover only 
the process of creation and decay of the pseudo-scalar rj meson, the so-called production 
current. However, in the complete process (1.2) the two nucleons can suffer initial state 
interaction (ISI) and final state interaction (FSI), before and after the rj creation, thus 
provoking distortions of the incoming and outgoing A'^A'^ waves. 

The ISI within a A^A^ pair before the rj creation is to be evaluated at relatively high 
energies, larger than the threshold of the rj meson production {Tkin ~ 1.3 GeV). Therefore, 
one can expect that the variation of ISI effects with the kinetic energy is small. Indeed, as 
shown in Ref . [35] , the effect of ISI can be factorized in the total cross section and it plays 
effectively the role of a reduction factor in each partial wave in the cross section. This 
reduction factor depends on the inelasticity and phase shifts of the partial waves at the 
considered energies. Near the threshold, the number of initial partial waves is strongly 
limited by the partial waves of the final states, and one can restrict oneself to ^Pq and ^Pi 
waves. Experimentally [36] it is found that at kinetic energies of the order of few GeV 
the phase shifts ^Pq and ^Pi are indeed almost energy independent and the reduction 
factor for each partial wave can be taken as constant. In our calculations we adopt for 
the reduction factor C, the expression from Ref. [35] 

= ^i{p) cos^ 5i{p) + ^iip)f , (4.1) 

where and 5i denote the inelasticity and phase shifts, respectively. We employ in our 
calculations C = 0.277 for the ^Pi wave (i = 1) and C = 0.243 for ^Pq (i = 2) [10]. 

FSI effects among the escaping A^A^ pair are accounted for within the Jost function 
formalism [37] which reproduces the singlet and triplet phase shifts at low energies, as 
appropriate for reactions near threshold. Details of calculations of FSI with the Jost 
function can be found in Ref. [38]. 
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B. T] Meson Production in NN Collisions 

Results for the energy dependence of the total cross section (j^'^-^^^v^ gq (3.30), are 
presented in Figs. 3 and 4, for proton-proton and proton-neutron reactions, respectively. 
The amplitude TArlr-vArAr?? represents a sum of nucleon, meson conversion and nucleon 
resonance currents, each of them being a sum over all the considered exchanged mesons 
(see Fig. 2). Using the parameters listed in Tabs. I and II the nucleon current contribution 
with parameters from the Bonn group [34] is found to be fairly small in the present 
calculations (see dot-dashed curves). 

Also the contribution from the conversion current (see doted curves) is not too large. 
The main contribution to the cross section comes from the nucleon resonance currents. 
Here it is worth mentioning that, in spite of the large number of the considered diagrams 
and the large number of the effective parameters, there is not too much freedom for fine- 
tuning of the cross section. As mentioned, most parameters are restricted by independent 
experiments and they cannot be varied in large intervals. We can slightly modify the 
less known parameters to achieve improvement of the overall description with data. In 
particular, in the present work we find a small contribution of the resonance £)i3(1520) at 
small energies, but a rather strong energy dependence due to tt and 77 exchange diagrams 
with increasing energy. Therefore, in order to reduce the influence of the 1^13(1520) at 
large excess energies, for these diagrams the cut-off parameters have been chosen smaller 
than for others (see Tab. II). 

Note that, even achieving a good description of the cross section in proton-proton re- 
actions, it is not a priory clear that the obtained set of parameters equally well describes 
also the proton-neutron reactions. The isospin dependence of the amplitude is determined 
by a subtle interplay of different diagrams with different exchange mesons (scalar, vector, 
isoscalar, isovector). Once the parameters for the pp reaction are fixed, the pn amplitude 
follows directly without further parameters (ISI and FSI are different for pp and pn sys- 
tems, but fixed independently) . Figure 4 demonstrates that the isospin dependence of the 
amplitude is correctly described with respect to available data. 

In Fig. 5 the angular distribution of 77 mesons in the center of mass is presented for two 
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excess energies, As^/^ = 16MeV (left panel) and As^/^ — 37 MeV (right panel). Here a 
few comments are in order. Generally, at threshold-near energies the angular distributions 
are rather flat. That means that the total cross section at the corresponding energy 
is 0"*°* oc Airda/dfl* and a quantitatively good description of the angular distribution 
da/dVt* provides, of course, also a good description of the total cross section. Therefore 
we compare our calculations with data at such an energy, where we achieve a better 
description of the total cross section (cf. Fig. 3). The other existing experimental data 
at similar energies (in particular, As^/^ = 15.5 MeV, [39]) differ from each other by an 
overall normalization factor, which can amount up to 40% at these energies, see Ref. [40]. 
Another important issue for the analysis of the angular distributions is the shape at 
forward and backward directions. The conversion current provides an upward (U) shape 
of the distributions, while the nucleon and nucleon-resonance currents generally give a 
downward (fl) curve at forward and backward directions. Since in our calculations the 
contribution of the conversion current is much smaller than the one from the nucleon 
resonance current, the resulting angular distribution, albeit being rather flat, has however 
an upward (U) trend. 

In Fig. 7 the invariant mass distribution of the two final protons is exhibited together 
with experimental data at an excess energy of Asi/2 = 15.5 MeV. Since our total cross 
section has been optimized to describe the data at higher energies, As^/^ > 16 MeV, see 
Fig. 3, the calculated curve in Fig. 7 has been multiplied by a factor of 1.4 obtained from 
an attempt to reconcile the data at As-*^/^ = 15.5 MeV and As^^'^ — 16 MeV (relevant 
comments and discussion about the normalization of data can be found in Refs. [39, 40]). 

The mild discrepancy of the angular distribution in Fig. 5 and the slight underestimate 
of the total cross section data at As^^^ ~ 10 MeV let us argue that further reaction details 
could be accommodated. For instance, the ss shake-off considered in [41] may be included. 
This, however, deserves a separate consideration. As interim summary we believe to have 
at our disposal an appropriate description of the rj production cross section. 

As mentioned in the introduction, in heavy-ion collisions a substantial part of di- 
electrons stem from Dalitz decays of rj mesons. In the same invariant mass region also 
A Dalitz and bremsstrahlung processes contribute to the total yield. It is, therefore. 
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important to have some reliable cross sections for elementary 77 production in pp and pn 
collisions as input for transport model simulations. Here one faces the problem whether 
the free cross sections (say, simply parameterizations of data) should be utilized or such 
ones where the FSI is "switched off", as in a dense hadronic environment the outgoing 
nucleons are not asymptotically free out-states. The same holds for the ISl. To have some 
indicator of the size of ISI and FSI effects we show in Fig. 6 the same as in Fig. 3 but with 
FSI switched off (dashed curve) or both ISI and FSI switched off (dotted curve, which 
may be called pure production cross section). A similar pattern holds for pn — > pnr] 
(not displayed). One observes indeed sizeable effects of ISI and FSI which point to the 
importance of a proper treatment of such effects in accurate many-body simulations. 



V. THE COMPLETE REACTION pp ^ ppr? ^ pp7e+e- 

In order to emphasize the dependence on the di-electron invariant mass we rewrite the 
integrated cross section (3.29) da/ds-y* — J dsj, {da / dsj,ds^*) in the form 

-'"|f'.r(v)P^^^^^X(v). (5.1) 



where we introduced dimensionlcss variables a, b and ^ 

da/ds^* may be considered as normalized contribution to the di-electron invariant mass 
spectrum mediated by the r] Dalitz decay. 

Since the total width of the 77 meson is fairly small, the parameter b in (5.2) provides 
a very sharp maximum of the integrand function at ^ = a as long as the parameter obeys 
a > 1 (or, equivalently, the di-electron invariant mass fulfills s* < rri^). This allows to pull 
out the smooth function cr^MN^NNri^^ — ^) from the integral and to calculate it at = m^, 
i.e. for the on-shell rj meson. The remaining part can be calculated analytically. However, 
for a < 1, when the di-electron mass is larger than the rj pole mass, the integrand does 
not exhibit anymore a resonant shape and the integral ought to be calculated numerically. 
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Equation (5.1) shows that the di-electron invariant mass distribution dajds.^* is pro- 
portional to the transition form factor F^^*{s^*) so that measurements of this distri- 
bution provide direct experimental information about the FF. It can be checked that 
up to a ~ 1.035 (corresponding to s^* ~ 0.29 GeF^) the smooth part of the integrand, 
^ArAr-»iVAr?7^^j^ pulled out from the integral at ^ = a obtaining 

X(V) - ^^^^ = (5.4) 

Equations (5.1 - 5.4) allow then to define an experimentally measurable ratio which is 
directly proportional to the form factor squared 



d(J I ds^* S^* /l S.y* ^rnin/ 



{dcT / ds^*)\^^^_^^^ Sj* ^jnin \ 1 ^7*/^^ / |-^?J77* (S7*,mm) P 

where s^*^rnin is the minimum value accessible experimentally (in the ideal case this is the 
kinematical limit Sj*^jnin — 4/Xg with electron mass /le)- At low enough values of s^*^rnin, 
the transition FF is close to its normalization point F^y*{0) = 1 and the ratio (5.5) is 
just the transition FF function of s^* . 

As a approaches unity keeping the maximum position still within the integration range, 
one can again withdraw from the integral the smooth function cr^^^^^n ^t — m^. 
However, now the function (^ — 1)^/^ can not be considered as smooth enough and must 
be kept under the integration. Nevertheless, even in this case the remaining integral can 
be computed analj^ically (see Appendix B) and one can still define a ratio analogous to 
(5.5) which allows for an experimental investigation of the FF near the free rj threshold 



'7* 



i^^). At > the integral X(s^*) does not have anymore a sharp maximum 
and it must be calculated numerically as mentioned above. 

In Fig. 8 the di-electron invariant mass distribution da/dslj^ is exhibited as a function 
of the invariant mass slj^ calculated from formula (5.1). The solid curve employs the 
transition FF F^^^. (s^.) calculated within the VMD model [5] 

with standard values for rUp and F^ [19]. (Note that in the kinematical region of interest 
the p contribution is sufficient.) The dashed curve in Fig. 8 is for the pure QED part, i.e. 
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without accounting for the strong form factor, i.e. with F^j^*{s^*) — 1. The experimental 
data in Fig. 8 are prehminary results from HADES [6] . It is clearly seen that without FF 
the calculated cross section rapidly drops as sl^i^ m^, while inclusion of the FF only 
mildly modifies the shape (factor 2 increase at m^* ~ 400 MeV High precision data are 

needed to arrive at a firm conclusion on the validity of the employed VMD FF. 

1/2 

With increasing s^* , the effect of the FF becomes more and more pronounced. This 
can be visualized if one calculates the e+e~ mass distribution beyond the free rj meson 
pole mass, i.e. at slf*^ > rrij^. Figure 9 exhibits the behavior of the mass distribution at 
large values of slj^ . The effect of the FF increases noticeably at large values of s^*. This 
is of interest since, as mentioned, the VMD calculations provides a good description of the 
on-shell t] meson FF in the range s^* < rn^. An investigation of the FF at s^. > can 
provide information of the relevance of the VMD model for off-shell hadrons at kinematical 
limit. In particular, when calculating Dalitz type processes with (off-shell) vector mesons 
created via bremsstrahlung off nucleon resonances, one usually estimates the couplings 
in the corresponding vertices from the radiative decay off an on-shell resonance with 
applying the VMD model. An additional justification of such a procedure can be obtained 
from investigating the transition FF beyond the on-mass shell limit. It should be noted, 
however, that severe background processes will make difficult an identification of the rj 
Dalitz yield in this kinematical region. 



VI. SUMMARY 



In summary we have analyzed the di-electron production from Dalitz decays of 77 
mesons produced in pp collisions at intermediate energies. The corresponding cross section 
has been calculated within an effective meson- nucleon approach with parameters adjusted 
to a large extent to describe the free vector meson production in A^A^ reactions near the 
threshold. We argue that by studying the invariant mass distribution of the final e+e" 
system in a large kinematical interval of di-electron masses, one can directly measure the 
77 meson transition form factor in, e.g., pp collisions. Such experiments are already 

performed and will be further improved at HADES. Our results may serve as prediction 
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for these forthcoming experiments. The uncertainties of such a procedure depend upon 
the scale of the background processes and is expected to be small if the interference is 
destructive [18]. Experimental information on transition form factors is useful for testing 
QCD predictions of hadronic quantities in the non-perturbative domain. 
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APPENDIX A: COUPLING CONSTANTS FOR SPIN-i RESONANCES 



The coupling constants of the spin-| resonance Di3(1520) can be estimated from the 
known decay widths. For instance, for the decay Di3(1520) Nrj one has 



Stt 



?77^ 4- 

V spins 



X Tr 



(Pat + mN)j5{PN* + rnN*)j5) 



(Al) 



where for the spin-| particles we use the Rarita-Schwinger spinors U^{Pn*-,s) for the spin 
s summation with the relation 



Y,W{Pn^,s)W{Pn*,s) = {PN-.+mM-^) 



X 



-9ny + 



7 7 ,17 ^N* 7 ^N*) , ^J^^N''^N* \ ^^2) 



+ 



3m jv* 



+ 



We use the short hand notation p — (p • 7). Equation (Al) provides slightly larger 
coupling constants than that used in Ref. [10] in which, besides the decay widths, the 
coupling constants have been adjusted to fit also other processes involving £)i3(1520). An 
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analog equation hold for Di3(1520) Ntt. In the present paper we take g^NN* — 1-55 
and Qr^NN* — 8.3, see Tab. II. 



APPENDIX B: INTEGRATION OVER THE VARIABLE s^y 

Since the parameter h in eq. (5.2) is very small in the whole kinematical region of s-y* 
the integrand has a sharp maximum around ^ = a. If a > Xmin = 1 then the maximum is 
located inside the integration interval of ^ and one can take advantage of the smoothness 
of the T] production cross section (j^^^^^v iq withdraw it from the integral. Then 
XC-af b/a 



MO = J dC- 



(Bl) 



= -C+ , f- In e + ^ In ((e - af + b') (2a' - 3a' + 2ab' - 3b' + l) 

a a{a^ + b^) I 2 ^ ^ ^ ^ 

+ ^arctan (^^^^ " + 3a'-a- 3ab' - b^ + 36^) . 

It is seen that in case when a ^ 1, 6 <^ 1 the term b^/a and terms containing logarithms 
in (Bl) can be disregarded: 

MO ^ ^^^arctan (^^) , (B2) 
(a - If 

lli^max) - lliirnin) ^ TT ^ • (B3) 

a 

This corresponds to the case when in eq. (5.2) one withdraws from the integral both the 
production cross section cr^^-^^^'n and the combination (^ — 1)^/^ at ^ = a. Note that 
eq. (Bl) is a good approximation even for a < 1, if a + 6 > 1, i.e. for a in the b vicinity 
of unity. 
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FIG. 1: Illustration of the choice of the independent variables for the process Ni + N2 — 
A^( + A''2 +7 + ^1 + ^2 (^1 = e"*", I2 = e~) within the duplication kinematics [25]. The invariant mass 
squared of the two final nucleons is denoted by S12, while the invariant mass of the subsystem 
7e"^e~ is s^. 




a) b) 

FIG. 2: Diagrams for the process NN — NN^e^e~ within an effective meson-nucleon theory, 
a) Dalitz decays of r] mesons from bremsstrahlung like diagrams. The intermediate baryon N* 
(triple line) can be either a nucleon or one of the nucleon resonances S'ii(1535), Pii(1440) or 
1)13(1520). Analog diagrams for the emission from Fermion line N2- b) Dalitz decay of r/ mesons 
from internal meson conversion. Exchange diagrams for identical nucleons in exit channel are 
not displayed. 
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FIG. 3: Total cross section for the reaction pp ppr] as a function of the excess energy 
As^/^ = ^/s — 2m — m^. The dot-dashed hne corresponds to the nucleon current contribu- 
tion (without resonances, Fig. 2a), while the dotted line is for the internal meson conversion 
diagram (Fig. 2b). The dashed curve exhibits the contribution of the S'ii(1535) resonance (nu- 
cleon resonance current in Fig. 2a). Contributions from Pii(1440) and Di3(1520) are smaller 
and not displayed. The solid curve is the total contribution with all interferences. Data are 
from Refs. [39, 40, 42-45]; error bars are suppressed. 
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FIG. 4: The same as in Fig. 3 but for the reaction pn pnr]. 
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FIG. 5: Angular distributions of rj mesons in the pp center of mass system for the excess energies 
As^/^ = 16 MeV (left panel) and 37 MeV (right panel). The dotted lines depict the contribution 
of the conversion currents, while the solid lines correspond to the total differential cross section, 
including also nucleon and resonance currents and interferences. Data are from Ref. [39]. 
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FIG. 6: Effect of ISI and FSI for the reaction pp ppij. The dashed curve depicts the cross 
section when switching off the FSI, while the dotted curve is without both FSI and ISI. The 
solid curve and data are as in Fig. 3. 
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FIG. 7: Distribution of the proton-proton invariant mass S12 for the reaction pp pprj at 
excess energy 

Asi/2 = 15.5 Mey. Since there is an overall normalization uncertainty of data 
at ^ 5 _ 16.0 Mel/ obtained by different groups our curve has been normalized to the 

data (see Ref. [40, 45], normalization factor « 1.4). Data are from Ref. [45]. 
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FIG. 8: Invariant mass distribution of di-electrons produced by Dalitz decay of r/ mesons in the 
reaction pp — > pp'je^e^ at initial kinetic energy Tfcj„ = 2.2 GeV . The sohd line is for results with 
the r/77* transition FF computed within the VMD model, while the dashed line corresponds to 
a pure QED calculation of the r/77* vertex, i.e., for {F^/'y-y* p = 1. The position of r] pole mass is 
exposed. Preliminary experimental data are from Ref. [6]. 
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FIG. 9: The same as in Fig. 8 but in the kinematical region beyond the "r/ threshold", i.e. 
for di-electron masses larger than the rj pole mass being accessible for intermediate off-shell r] 
mesons. 



30 



